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ABSTRACT 

We compare the galaxy population of a smoothed particle hydrodynamics (SPH) simu- 
fS ^ lation to those predicted by the GallCS N-body + semi-analytic model and a stripped 

down version of GallCS that omits supernova and AGN feedback. The SPH simula- 
tion and the no-feedback GallCS model make similar predictions for the baryonic mass 
functions of galaxies and for the dependence of these mass functions on environment 
and redshift. The two methods also make similar predictions for the galaxy content 
of dark matter haloes as a function of halo mass and for the gas accretion history of 
galaxies. There is a fairly good correspondence between the "cold" and "hot" accretion 
Jv>( \ modes of the SPH simulation and the rapid and slow cooling regimes of the GallCS 

■ calculation. Both the SPH and no-feedback GallCS models predict a bimodal galaxy 

' population at z = 0. The "red" sequence of gas poor, old galaxies is populated mainly 

by satellite systems, which are starved of fresh gas after they begin orbiting in larger 
haloes, while, contrary to observations, the central galaxies of massive haloes lie on 
the "blue" star-forming sequence as a result of continuing hot gas accretion at late 
times. Furthermore, both models overpredict the observed baryonic mass function, 
especially at the high mass end. In the full GallCS model, supernova-driven outflows 
reduce the masses of low and intermediate mass galaxies by about a factor of two. 
AGN feedback suppresses gas cooling in large haloes, producing a sharp cut-off in the 
baryonic mass function and moving the central galaxies of these massive haloes to the 
red sequence. Our results imply that the observational failings of the SPH simulation 
and the no-feedback GallCS model are a consequence of missing input physics rather 
than computational inaccuracies, that truncating gas accretion by satellite galaxies 
automatically produces a bimodal galaxy distribution with a red sequence, but that 
explaining the red colours of the most massive galaxies requires a mechanism like AGN 
feedback that suppresses the accretion onto central galaxies in large haloes. 

Key words: cooling flows — galaxies: evolution — galaxies: formation — galaxies: 
haloes — galaxies: ISM 



1 INTRODUCTION 

The classification of galaxies into early-type and late-type 
goes back to the earliest studies iHubblelll926l : iHumaso-r] 



Il936l) . but the clear demonstration that there is no smooth 
continuity between the two types wa s only possible with 
the advent of large statistical surveys dStrateva et aljfeuOll : 



Cattaneo et al. 



Kauffmann ct al. 


l200l 


2004: iHoee et al. 


120041: 



Bald rv et all [2004; Ba logh et alJ 



a comprehensive introduction to the bimodality problem). 
Galaxy colour-magnitude diagrams from the Sloan Digi- 
tal Sky Survey (SDSS) show a clear distinction between 
a tight red sequence with no recent star formation and a 
broader blue sequence, where star formation is still going 
on. The blue sequen ce is truncated n ear the characteristic 
luminosity L* of the ISchechterl ]l976lj luminosity function, 
while the red sequence extends to luminosities several times 
higher. Elliptical galaxies concentrate in the bright part of 
the red sequence (L > L»). Galaxies in crowded environ- 
ments are normally red, while blue galaxies reside main ly in 
the field (iKauffmann et al J 1200 A iBlanton et alJ 120051: and 
numerous references therein). Deep surveys like COMBO- 
17 and DEEP- 2 show that bimodality is already well estab- 
lished at z ~ 1 jBell et all2004l ; lFaber et all200ai . and even 
at 2 > 3 there appears to be a distinction between actively 
star-forming galaxies and pa ssive systems dominated b y red- 
der stellar populations (e.g.. Ivan Dokkum etaiT 2006 ) . Here 
we investigate the physical origin of this bimodality in the 
context of modern galaxy formation theory, comparing the 
results of two different modelling methods. 

Cosmological hydrodynamic simulations (see 
iFrenk et al.lll999l for a review), which integrate the equa- 
tions of motions for the baryons alongside those for the dark 
matter, and hybrid simulations, which combine an N-body 
treatment of the d ark matter with sem i -analytic modelling 
of the baryons (e .g. |Yo^hida^et_al: |200i|Hattoii^^lJ 

ICroton et alJ 



Hellv et al]l2003l IPeLucia et alJ 



2004: 



2003; 



2005; 



Cattaneo et al . 2006), are the main methods for exploring 



the physics of galaxy formation in a cosmological scenario 
while retaining the spatial information on the galaxy distri- 
bution. The fully numerical approach treats the dynamics 
and the gas physics more self-consistently because accretion 
and mergers come out naturally from the physics of the 
simulations. In the hybrid approach, the same phenomena 
are modelled through assumptions, such as the hypothesis 
that the gas starts cooling from a spherical distribution at 
the virial temperature of the halo or the dynamical friction 
formula to compute merging rates. However, hydrodynamic 
simulations are expensive in CPU time. That limits the 
resolution mass and the size of the computational box, the 
latter resulting in poor sampling of rare objects. It also 
means that, for the same CPU time, hybrid methods allow 
one to explore a wider set of physical assumptions, such 
as different feedback scenarios. In this paper, we compare 
results fro m a smoothed p article hydrodynamics (SPH) 
simulation llKatz et al.lll996T) to those of the GallCS hybrid 
mode l (Galaxies In Cosmological Simulations: lHatton et all 
120031) . 

This comparison tests the consistency between the gas 
dynamical treatment of the SPH simulation and the ide- 
alised model incorporated in GallCS, since we run Gal- 
ICS with merger trees extracted from the SPH simula- 
tion i tself. Two groups , centred in D urham jBejrscj ietalJ 
120011: iHellv et alJl2003D and Munich (lYoshida et alJl2002t) . 
have already performed similar checks, verifying that the 
two semi-analytic mod els investigated llCole et alJ l2000t 
IKauffmann et al] Il99^ predict galaxy masses and lumi- 
nosity functions in reasonable agreement wit h cosm ologi- 
cal SPH simulations run with the HYDRA dPearce et al] 



1200 ill and GADGET dSpringel et al.l l200ll) codes, when 
the physical assumptions are similar. One motivation for 
our study was simply to conduct the same test for Gal- 
lCS. However, we also focus on issues related to the 
galaxy bimodality, such as the distinction between cen- 
tral and satellite galaxies and the possible connection 
to t he distinction between "cold" and "hot" gas accre- 
tion llKatz et al.ll2003l:lBirnboim fc DekejEoolt iKeres etafl 
120051 : |p'ekel^^BirriboiirI]2^0^) ~ To illuminate these issues, 
we consider two versions of the GallCS model: a "stripped- 
down" version with no feedback, designed to approximately 
match the p hysics of the SPH sim ulation, and t he full version 
descri b ed bvlHat ton et all ]2003|) a nd used in|Blaizo^^fl] 
]2004h . lLanzoni et al] ]2005f) and ICattaneo et al.l ]2005l) 
which contains supernova and AGN feedback prescriptions 
designed to reproduce the blue luminosity function of galax- 
ies. 1 We show that the SPH simulation and the no- feedback 
GallCS model produce a red sequence of satellite galaxies in 
massive haloes, but predict, contrary to observations, that 
the central galaxies of massive haloes lie on the blue se- 
quence. The strong supernova feedback incorporated in the 
full GallCS model reduces galaxy baryonic masses to a level 
consistent with observations. The combination of supernova 
and (more importantly) AGN feedback shuts off gas accre- 
tion in massive central galaxies, turning them red. 

Section 2 describes the astrophysical and computational 
assumptions used in the SPH simulation together with the 
algorithms used to identify galaxies and haloes, to attribute 
galaxies to haloes, and to construct merger trees by link- 
ing haloes identified at different time-steps. Section 3 sum- 
marises the GallCS implementation used here, in particu- 
lar its assumptions about gas accretion, star formation and 
feedback. In Section 4 we begin our comparison of the two 
methods with a global view of the model galaxy popula- 
tions. We investigate the mass function, its dependence on 
redshift and environment, and the relation between galaxies 
and dark matter haloes. We then move into more specific 
aspects such as the growth mechanism (cold flows versus 
hot flows; Section 5), the gas content, and the presence (or 
absence) of a bimodality in the distributions of gas content 
and star formation time-scales (Section 6). In Section 7, we 
discuss how this work contributes to our understanding of 
galaxy formation and of the processes that are at the basis 
of the observed galaxy bimodality. Section 8 summarises our 
results. 



1 As we discuss later, a newer version of GallCS, with an im- 
prov ed treatment of AGN and supernova feedback, is now avail- 
able JCattaneo et all2006t) . The new version can fit the joint dis- 
tribution of galaxy magnitudes and colours at 2 ~ and the 
luminosity function of Lyman-break galaxies at 2 ~ 3 to unprece- 
dented accuracy. However, we carried out all of our analyses prior 
to the development of this new version, which was partly moti- 
vated by ideas discussed here. The old version is adequate for 
the purposes of this paper. We also wanted to provide tests of 
the GallCS version used for most of its results that have been 
published to date. 
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2 THE SPH SIMULATION 

2.1 Cosmological model and simulation 
parameters 

The cosmology used for this study is a flat ACDM uni- 
verse (inflationary cold dark matter with a cosmological con- 
stant) with h = ,Ho/(100 kms* 1 Mpc" 1 ) = 0.65, Q, m = 0.4, 
Q.A = 0.6, tt b = 0.02/i~ 2 = 0.0473, inflationary spectral in- 
dex n = 0.93, and power spectrum normalisation cr$ = 0.8. 
We have used this combination of parameters in a series 
of simulation s over the course of a number of years (be- 
ginning with iDave et al 1 ll999l) . systematically varying nu- 
merical resolution and simulation volume. The values of 
h, fit, n, and as are similar to those inferred from joint 
analyses of cosmic microwave background (CMB) measure- 
ments from the WMAP satellite and galaxy power spec- 
trum measurement s from the 2dF Galaxy Redshift Survey 
and theSDSS fe.g JSnergel et al.teOOfllT^rnark et al.l20oi 
ISanchez et al.l200a) . These analyses generally favour a lower 
Qm value, in the range ~ 0.23 — 0.3, but the qualitative fea- 
tures of galaxy formation are not sensitive to the precise 
value of tt m (a point we have checked by running a compa- 
rable simulation with Q m = 0.3). Since we adopt identical 
parameters for the SPH and GallCS modelling, this simu- 
lation should be entirely adequate for our purposes in this 
paper, despite its slightly outdated parameter values. 

The simulation volume is a comoving periodic cube 
22.22ft" 1 = 34.19 Mpc on a side, modelled using 128 3 dark 
matter particles and 128 3 gas particles. Gravitational forces 
are softened with a cubic spline kernel of comoving ra- 
dius 7.7 kpc, approximately equivalent to a Plummer force 
softening of e gra v = 5.4 kpc, while hydrodynamic pressure 
forces are calculated over 32 neighbours. The baryonic mass 
threshold for resolved galaxies is 6.8 x 1O 9 M0 , the mass of 64 
gas particles, and there are 1120 galaxies in the box above 
this thresho ld at z = 0. Furthe r details of the simulation can 
be found in lKeres et all l)2005h . 



2.2 The SPH code 

The simulatio n was performed with t h e parallel version 
of TreeSPH iHernauist fc Katzl Il989t iKatz et al.l Il99fit 
IDave et alJ Il997l) . This code combines smoothed parti- 
cle hy drodynamics ( SPH: lLucvlll977l : iGingold fc Monaghanl 
ll977T with a hiera rchical tree algorithm for computin g 
gravitational forces fearnes fc Hudll986l : lHernauisdll987t) . 
TreeSPH is a completely Lagrangian code, adaptive both in 
space and in time. It contains three kinds of particles, rep- 
resenting dark matter, gas and stars. The dark matter and 
the stars are only subject to gravity, while the gas is also 
subject to pressure gradients and shocks. The gas experi- 
ences adiabatic heating and cooling, shock heating, inverse 
Compton cooling off the microwave background and radia- 
tive cooling via free-free emission, collisional ionisation, col- 
lisional recombination, and collisionally excited line cooling. 
The cooling rate is calculated for primordial chemical abun- 
dances and, since there is only atomic cooling, the gas cannot 
cool radiatively below T ~ 10 4 K. A uniform photoionising 
UV background heats low temperature gas and suppresses 
cooli ng processes involving n eutral atoms at low gas densi- 
ties teaardt fc Madadll996h . 



The simulation also includes a simple treatment of star 
formation and its associated supernova feedback. Gas above 
a threshold hydrogen number density nu = 0.1cm -3 turns 
into stars on a time-scale set by the dynamical time or the 
cooling time, whichever is longer. Additional star formation 
conditions are that the gas is Jeans unstable, part of a con- 
verging flow and above the virial overdensity, but gas that 
satisfies the density criterion usually also satisfies the other 
three. Gas reaches this high density only after cooling to 
T ~ 10 4 K. The subsequent molecular and metal-line cool- 
ing to lower temperatures is implicitly included as part of 
the star formation process. This star formation prescription 
leads to a rel ation with the gas surface dens ity similar to a 
Sc hmidt law |Schmiddh95 9|:|Kennicutlll998l see discussions 
bv lKatz et al.lll99rj and lStinson et al.ll2006l) . 

Stars more mas s ive th an 8Mq explode as supernovae. 
For a lMiller fc S calo ( 1979) initial mass function that means 
7.35 x 10 supernovae per solar mass of formed stars. Each 
supernova releases 10 erg, gradually distributed as heat to 
the gas particles nearby with an exponential time decay of 
2 x 10 7 yr. The surrounding medium is usually dense, so 
the deposited energy is typically radiated away before it can 
drive a galactic scale wind. For this reason, we consider this a 
"minimal" feedback algorithm, even though the full amount 
of expected supernova energy is incorporated. 

2.3 Identification of galaxies 

Cosmological simulations that incorporate cooling and star 
formation produce dense groups of baryonic p articles with 
the sizes and ma sses of observed galaxies dKatzl 119921 : 
lEvrard et all Il994h . We identify these aggregations with 
the group finding a lgorithm 'Spline Ker nel I nterpolative 
DENMAX' (SKID) 2 iGelb fc BertschingeJl994:lKatz et all 
1996). This algorithm involves five basic steps: (1) determin- 
ing the smoothed baryonic density field; (2) moving baryonic 
particles towards higher density along the initial gradient of 
the baryonic density field; (3) defining the initial group as a 
set of particles that aggregate at a particular density peak; 
(4) linking together initial groups that are very close to- 
gether; (5) removing group particles with positive binding 
energy relative to the group's centre of mass. We apply SKID 
to the population of all star particles and to the gas parti- 
cles that have temperatures T < 3 x 10 4 K and overdensities 
Pgas/pgas > 10 3 . Henceforth, we refer to the aggregations of 
stars and cold gas that SKID identifies simply as "galaxies." 
Tests on simulations with varying mass resolution show that 
the simulated galaxy population becomes substantially in- 
complete below a baryonic mass corresponding to ~ 64mspH 
but i s fairly robust above this limit (see, e.g.. iMurali et all 
,2002). We therefore adopt 64m S p H (6.8 x 1O 9 M ) as our 
resolution threshold and ignore lower mass galaxies in our 
analysis. Because of our high overdensity threshold for star 
formation, essentially all star formation in the simulation 
takes place in galaxies, though some of these are below the 
resolution limit, and some stars are tidally stripped from 
galaxies during dynamical interactions. 

The minimum galaxy mass of 64mspH corresponds to 

2 We use the implementation of J. Stadel and T. Quinn, available 
at http : //www-hpcc . astro .Washington. edu/tools/ skid.html 
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a minimum host halo mass of ~ (f2 m /f2;,)64mspH, since the 
fraction of cold gas in a halo never substantially exceeds 
the universal baryon fraction. The corresponding minimum 
virial temperature for our simulation is ~ 80,000(1 + z) K 
for the simulation considered in this article, where the (1 + z) 
factor arises from the increasing physical density at higher 
z for fixed virial overdensity. The redshift dependence of the 
virial overdensity makes the '80, 000 K' factor slightly lower 
at high redshift (~ 75, 000 K) and slightly higher at low red- 
shift (88, 000 K at z = 0). Any resolved galaxy in the simu- 
lation resides in a dark matter halo with virial temperature 
higher than this minimum temperature. 



3 THE GALICS HYBRID MODEL 

The semi-analytic model that we use in t his paper is exactly 
the one described in lHatton et, alJ (|2003ft . In this section, we 
therefore only briefly summarise those of its features which 
are relevant to the present study. 

3.1 Dark matter haloes and their merger trees 

We identi fy and charac t erise dark matter haloes the same 
way as in (Hatton et alJ J2003[). We use a friends-of-friends 
(FOF) algorithm l)DavisetHflTl98fift to detect groups of par- 
ticles with overdensity ~ 200 times the mean density (using 
a fixed linking-length parameter b — 0.2). For each group 
with more than 20 particles, we compute the total kinetic 
and potential energies, and only keep in the halo catalogue 
the groups that are bound. The halo mass we refer to in the 
rest of the paper is simply N p x m p , where N p is the number 
of particles in the group, and m p the mass of a dark matter 
particle. The minimum halo mass is thus ~ 1.8 x 10 10 M Q . 
We set the centres of groups to be at the positions of their 
centres of mass. 

The information about an individual halo measured 
from the SPH outputs and passed to the semi-analytic model 
is contained in three parameters: the virial mass, the spin 
parameter, and the virial radius. The virial mass and the 
virial radius are measured in two steps. First, we compute 
the inertia tensor of each halo. From this we determine the 
three main axes of an ellipsoid that fits the mass distribution 
and is centred on the centre of mass. Then we shrink this 
ellipsoid until the halo particles that it contains satisfy the 
virial theorem. The mass of the particles in the ellipsoid at 
this point is the virial mass. The virial radius is the radius 
of a sphere with the same volume as the virial ellipsoid. 

We compute merger trees by linking haloes identified in 
each SPH snapshot with their progenitors in the previous 
one. All predecessors from which a halo has inherited one or 
more particles are counted as progenitors. 

3.2 The cooling scheme 

Newly identified haloes receive a gas mass determined by 
the universal baryonic fraction fi(,/fi m . As in other semi- 
analytic models, all baryons start in a zero-metallicity hot 
phase, shock-heated to the virial temperature. We assume 
that the the hot gas density profile is described by a singular 
isothermal sphere truncated outwards at the virial radius 
and inwards at a core radius of 0.1 kpc. 



The cooling time t C ooi(r) is calculated from the hot 
gas density di stribution with the metal-de pendent cooling 
function from ISutherland fc Dopital <ll993ft . There are two 
differences with the SPH calculation here. First, the semi- 
analytic calculation does not include a photoionising back- 
ground, although photoionisation effects on galaxy forma- 
tion ar e small in the halo mass range resolved by the simu- 
lation jQuinn et al.ll996tlThoul fc Weinbersll996ft . Second, 
GallCS does take into account the enhancement of cool- 
ing by metals, contrary to the SPH code, which assumes 
zero-metallicity cooling. However, in the no-feedback Gal- 
lCS model, metals are never expelled from galaxies, and the 
intra-cluster medium always remains of pristine composi- 
tion. Hence, the stripped down version of GallCS reproduces 
a cooling scheme that is in practice very similar and hence 
directly comparable to that of the SPH simulation. 

In small haloes, the cooling time is very short and all 
the gas gets cold almost immediately. Without thermal pres- 
sure support, the gas collapses to the centre of the dark mat- 
ter halo in free fall until the centrifugal force determined by 
the conservation of angular momentum balances the gravita- 
tional force and the gas settles into a disc. In this regime, the 
gas accreted by the disc in the time At between two snap- 
shots is the gas that is found within a radius Hnfaii = At/v c 
from the halo centre (here v c is the halo's circular velocity). 
In larger haloes the cooling time is longer, and gas slowly 
flows to the centre as it cools. In this case it is the cooling 
time that determines the rate at which the disc can grow. 
The gas that can cool in the time At between two snap- 
shots is the gas within a sphere of radius r coo i such that 
icooi(T'cooi) = At. The At interval between redshift outputs 
is typically ~ 0.3 Gyr (~ 0.1-ff -1 (z)) at high redshift and 
~ 0. 5 Gyr (~ QmH^i z) ) at l ow redshift. 

iBirnboim fc Dekell 12003ft . iKatz et alJ <2003ft . and 
iKeres et^dTl2005ft emphasise the distinction between "hot" 
gas accretion, in which shocks heat gas to the halo virial tem- 
perature near the virial radius, and "cold" accretion in which 
cold, unshocked gas streams penetrate far inside the virial 
radius. To examine this distinction in GallCS, we identify 
the cooling-limited and infall-limited regimes with the "hot" 
and "cold" accretion modes, respectively (see Section 5 be- 
low). The boundary between these regimes depends on At, 
and hence on our particular choice of the redshift outputs 
used for constructing halo merger trees. The physical identi- 
fication with_the_two accre tion modes is therefore on ly qual- 
itative. IKeres et alJ i2005ft and ICroton et alJ <l2005ft discuss 
the relation between the accretion modes and the cooling 
criteria of semi-analytic models in greater detail. 

3.3 Star formation 

Star formation is activated when the gas surface density is 
Egas > 20m p cm~ 2 , where m p is the proton mass. Disc gas 
forms stars at a rate M* = M co id/i*, where M co id is the 
mass of the gas availab le to form stars and t* = 50 tdyn 
iGuiderdoni et al]ll998l : the dynamical time tdyn is the time 
in which a star completes half a circular orbit at the disc 
half mass radius). The disc is assumed to have an exponen- 
tial profile and its radius is determined from conservation of 
angular momentum. 

Mergers and disc instabilities form bulges by transfer- 
ring gas and stars from the disc to the spheroidal compo- 
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nent faatton et alj l2003). In this paper, we are not directly 
concerned with morphologies. However, these morphological 
transformations can temporarily accelerate the conversion of 
gas into stars, and the bulge mass determines the effective- 
ness of AGN feedback as described below. 

3.4 Feedback 

The formation of young stars is rapidly followed by super- 
nova explosions. The energy released by supernovae is used 
to reheat the cold gas and remove it from the galaxy. The 
vers ion of the GallCS m odel used to produce this paper fol- 
lows lHatton et all ((2003) and uses a feedback prescription in 
which supernova outflows entrain more gas in massive galax- 
ies owing to the larger porosity of the interstellar medium 
( Silk 2001, 2003). This model leads to the very simple result 
that the outflow rate is comparable to the star formation 
rate, independently of the depth of the potential well. 

Supernova feedback is not enough to compensate the 
overcooling of gas in massive haloes, which would produce 
too many massive galaxies. Following ideas that have been 
around in the recent literature, we postulate the existence of 
a second feedback mechanism, associated with the growth of 
supermassive massive black holes in early-type galaxies. The 
simplest way to model black hole feedback is to assume an 
energy threshold. When the total energy injected by active 
galactic nuclei (AGN) in the intergalactic medium (IGM) of 
a group or cluster of galaxies exceeds the threshold, the hot 
gas is assumed to have acquired so much entropy that it is no 
longer able to cool. As black hole accretion releases an energy 
proportional to the mass accretion rate, and as the black hole 
mass is proportional to that of the host bulge, this criterion 
translates into a requirement on th e total mass of bulge stars 
in the group or cluster. Following lHatton et alJ i2003t) . we 
assume that gas stops cooling when £Mb u i gc > lO n M0, 
where the sum is over all galaxies in the halo. 



4 GLOBAL VIEW OF THE GALAXY 
POPULATION 

We consider three models: the SPH model, the GallCS 
model without any type of feedback, and the GallCS model 
with the feedback recipes described in Section 3.4. Because 
the feedback incorporated in the SPH simulation does not 
generally drive galactic winds or suppress gas cooling in mas- 
sive haloes, it is the no-feedback GallCS model that is most 
comparable to the SPH simulation in its physical assump- 
tions. 

Figure 1 compares the galaxy baryonic mass functions 
at z — in these three models. (Here baryonic mass in- 
cludes stars and the cold interstellar medium, but it does 
not include hot gas in the galaxy or group halo.) It also 
shows the separate contribution to the mass function of 
"field" galaxies (M ha i < 3 x 10 12 M Q ), " group" galaxies 
(3 x 10 12 M e < M ha i < 10 14 M Q ), and "cluster" galax- 
ies (M ha io > 10 14 M Q ). In fact, the 34.19 Mpc simulation 
cube contains only a single "cluster" mass halo (M = 
3 x 1O 14 M0), and the most massive "group" haloes are 
~ 3 x 10 13 M Q . 

The SPH and no-feedback GallCS predictions agree re- 
markably well, both for the total baryonic mass function 
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Figure 1. The galaxy mass functions computed with three dif- 
ferent models (the SPH model, the no-feedback GallCS model, 
and the GallCS model with the feedback reci pes of Sect i on 3.4 ) 
compared with the observational estimate of iBell et all J2003I) . 
In each of the three diagrams, the observational mass function is 
shown by the points and the theoretical one by the red line. The 
latter is decomposed into the mass functions of "field" galaxies 
(black solid line), "group" galaxies (dashed line) and "cluster" 
galaxies (dotted line), as described in the text. 



and for the separate mass functions in the three halo mass 
regimes. Feedback in the full GallCS model reduces galaxy 
masses, dropping the low mass end of the mass function and, 
most strikingly, producing sharp truncation of the high ends 
of the mass functions in the group and cluster haloes. 

Points show the observa tional estimate o f the total 
baryonic mass function from IBell et all ((2003) , based on 
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Figure 2. The galaxy content of dark matter haloes in the SPH simulation and in GallCS. The blue points correspond to galaxies that 
are at the centre of their dark matter halo, while the red points are satellite galaxies. In the SPH simulation, we always define the most 
massive galaxy to be "central." In GallCS, a halo may have no central galaxy if it has had a recent merger, because dynamical friction 
takes time to make galaxies sink to the new centre of mass. 



2MASS and SDSS data. Because of the limited simulation 
volume, we would not expect perfect agreement with the 
data even for a model with exactly correct input physics. 
Nonetheless, the SPH and no-feedback GallCS predictions 
clearly exhibit the well known tendency of models with effi- 
cient gas cooling and minimal feedback to overproduce the 
baryonic mass function, a dis crepancy that rem ains with 
larger simulation volumes fsee iMurali et all 2002). The full 
GallCS model, by contrast, appears to underpredict the high 
mass end of the mass function. We believe that this is a 
real discrepancy, arising because the lHatton et al] (I2003T) 



feedback prescriptions are tuned to match the _B-band lu- 
minosity function, and they produce massive galaxies that 
are too blue. The ne wer version of GallCS described by 
ICattaneo et al.l l)200fit) shuts off gas accretion onto central 
galaxies of all high mass haloes and thereby achieves a simul- 
taneous match to the u- and r-band luminosity functions. 

Figure 2 portrays the galaxy populations of all haloes in 
each of the three models, at z = 3 and 2 = 0. The "central" 
galaxy of each halo is marked by a blue point, while satellites 
in the same halo are marked by red points directly under- 
neath. In GallCS, as in other semi-analytic models, central 



Accretion, feedback and galaxy bimodality 7 




vn 11 12 13 14 15 W 11 12 13 14 15 

Log(M holo /M ) Log(M holo /M ) 

Figure 3. The fraction of the total baryonic content of a halo in stars (dashed lines) and in cold gas plus stars (solid lines) as a function 
of halo mass. The results of the SPH simulations are printed in black. The results of the GallCS model without feedback are printed in 
red. The results of the standard GallCS model are printed in blue. 
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Figure 4. The mass function of galaxies in the SPH simulation (points) and in the GallCS model (lines) at four different redshifts. In 
each panel, the dashed line is GallCS without supernova or AGN feedback, the solid line is the standard GallCS model, and the dotted 
line is the standard GallCS model when all masses are multiplied by a factor of 2.1. The mass function is given as the number of galaxies 
per logarithmic bin of 0.25 dex. In contrast to Fig. 1, we have not divided the number of galaxies by the volume of the computational 
box, so that the numbers convey an idea of the statistics in each bin. At Af ga i > 3 X 10 11 A/q, the number of galaxies in the box is < 10, 
and statistical effects (quantified by the Poisson error bars) are no longer negligible. 
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Figure 5. Mean halo occupations measured for the SPH galaxies (asterisks) and the full GallCS model galaxies (diamonds), at z = 0. In 
the left (resp. right) panels, galaxies where selected in baryonic mass so as to yield a number density of O.lfe 3 Mpc -3 (resp. 0.02h 3 Mpc -3 ). 
From top to bottom, we show the mean occupation functions for the full, old, and young populations. 



galaxies are treated specially — any cooling from the halo 
joins the central galaxy, not the satellites. The SPH sim- 
ulation does not impose this physics a priori, but central 
galaxies nonetheless emerge as a distinct class. The galaxy 
closest to the halo centre-of-mass is almost always the most 
massive galaxy in the halo and it us ually has the oldest 
stellar population feerlind et al.ll2003l) . Since the centre-of- 
mass is affected by asymmetries in the halo structure, the 
closest galaxy is occasionally a lower mass "satellite," with 
the most massive galaxy nearby. For convenience, we have 
simply defined the most massive galaxy of each SPH halo to 
be the central galaxy. On closer inspection, this identifica- 
tion appears reasonable in virtually every case. While every 
SPH halo thus has a central galaxy by definition, a GallCS 
halo may not, if the halo has experienced a recent merger 
and there has not been sufficient time for dynamical friction 
to drag the most massive galaxy to the bottom of the halo 
potential well. 

In all three models, the central galaxy baryonic mass 
is tightly correlated with halo mass, while satellite galaxy 
masses are broadly scattered with little correlation. Since 
most gas cooling in the SPH simulation and all gas cooling in 
the GallCS model occurs on the central galaxy, its growth is 
tightly coupled to the growth of the host halo. Most satellite 



galaxies, on the other hand, grew in smaller haloes, and 
they stopped growing after merging into a larger halo and 
becoming satellites. 

At A/ ha i < 3 x I0 12 M Q , both in the SPH model and 
in the GallCS model without feedback, most of the blue 
points are aligned very close to the M ga i = (fib/f2m)Mhaio 
black diagonal line, at least at z = 0. That means that most 
of the halo baryons are in the central galaxy. At Mhaio > 
3 x IO 12 M0, this is less and less true for two reasons: a) 
in dense environments the cumulative baryonic content of 
satellite galaxies is more important, and b) as the cooling 
time of the hot gas becomes longer, a larger fraction of the 
baryons remain in the hot phase (Fig. 3). For these reasons, 
we take Mhaio = 3 x I0 12 Mq as the border that separates 
field galaxies from group galaxies. 

The addition of supernova and black hole feedback has 
two important effects. At low halo masses, ejection of gas by 
supernova feedback shifts the locus of central galaxy points 
down by about a factor of two, with concomitant reductions 
in satellite galaxy masses. The factor of two emerges because 
the ejected mass is, by construction, about equal to the mass 
that forms stars. This suppression is evident both at z = 3 
and at z — 0. At z = 0, black hole feedback has a much 
more drastic effect on massive galaxies, introducing a sharp 
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cut-off at a total baryonic mass ~ 3 x 10 1 Mq. The most 
massive galaxies in the no-feedback model are up to a factor 
of ten larger. 

Figure 3 shows the mean mass of stars (dashed lines) 
and stars plus cold gas (solid lines) in bins of halo mass, 
now including central and satellite galaxy contributions. The 
SPH and no-feedback GallCS results agree fairly well at 
both z = 3 and z — 0, with GallCS predicting somewhat 
higher masses in the most massive z = haloes. The lower 
SPH stellar masses in low mass haloes are probably an arte- 
fact of limited numerical resolution, which leads to underes- 
timated star formation rates and stellar masses in galaxies 
near the resolution threshold. We again see the factor of 
two reduction in baryonic masses due to feedback in the full 
GallCS model, which grows to a larger factor at high halo 
masses. Stellar masses are suppressed by about the same fac- 
tor. Since the initial total gas content of each GallCS halo 
is (f2i,/f2 m )Mhaio by assumption, the baryons that are not 
in the cold gas plus stars phase are necessarily in the hot 
phase. 

Figure 4 compares the galaxy baryonic mass functions 
at z = 3, 1, 0.5, and 0. The generally good agreement be- 
tween SPH and no-feedback GallCS seen in Figure 1 extends 
to higher redshifts. However, galaxy growth in no-feedback 
GallCS lags slightly (an effect also seen in Fig. 2), so its mass 
function begins somewhat lower at z — 3 and catches up at 
low redshift. This difference may reflect the hig h efficiency 
of filamentary accretion fsee iKeres et alJl2005l) relative to 
the spherical accretion assumed in GallCS. The dotted lines 
in Figure 4 show the full GallCS mass function (solid line) 
after all galaxy masses are multiplied by a factor of 2.1. This 
shift brings it into good agreement with the SPH and no- 
feedback GallCS mass functions up to M ~ 3 x 10 M©, 
beyond which one sees the sharper cut-off induced by AGN 
feedback. 

Figure |S] compares the halo occupation statistics of 
galaxies in the SPH simulation (asterisks) and the full Gal- 
lCS model (diamonds), at z = 0. In each model, we select 
galaxies above a baryonic mass threshold that yields a mean 
space density of 0.1/i 3 Mpc" 3 (left panels) or 0.02/i 3 Mpc -3 
(right panels). We have divided the galaxies in each panel 
into an old and a young population based on the median 
age of their stellar population. The old/young threshold was 
chosen in such way that the two subsets contain the same 
number of galaxies. Top, middle, and bottom panels show 
the mean occupation function (the average number of galax- 
ies as a function of halo mass) for the full, old, and young 
populations, respectively. The SPH and GallCS predictions 
agree impressively well, including the separate halo occu- 
pation statistics of old and yo ung populations, wh ich are 
quite different from each other. iBerlind et alJ (I2003T) found 
a similar level of agreement between a larger volume, lower 
resolution SPH si mulation an d the s emi-analytic galaxy for- 
mation model of ICole et all (l200Cl) . Although we do not 
show it here, the SPH and GallCS models exhibit compa- 
rably good agreement for the mean pair counts and triple 
counts ((N(N-l)) and (N(N - l)(N - 2)}) as a function of 
halo mass. The agreement of halo occupation statistics es- 
sentially guarantees that the SPH method and GallCS will 
make similar predictions for large scale clustering, even in 
larger volumes that are more statistically representative, and 
that this agreement will extend to the dependence of clus- 
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Figure 6. Growth history of the central galaxy of the third most 
massive halo (A/halo ~ 3 X 10 13 Mq) in the SPH simulation (solid 
lines) and in the GallCS model without feedback (dashed lines). 
The black lines show the total galaxy mass (stars plus cold galac- 
tic gas), while the red lines show mass of the gas that has cooled 
onto the galaxy since the beginning of the simulation. Galaxy 
mergers are the reason why the black lines are above the red 
lines. The blue lines show the cold galactic gas content in the 
SPH simulation and in the GallCS model without feedback. 

tering on galaxy age and to the 3-point correlation statistics. 
On scales comparable to or smaller than the virial diame- 
ters of the largest haloes, the predicted spatial distribution 
of galaxies within haloes matters as well as the number of 
gal axies per halo . In in vestigating this issue, we found that 
the lHatton et all i2003T) prescription for computing satellite 
galaxy positions predicts spatial distributions in large haloes 
that are to o diffuse, artificially suppressing small scale clus- 
tering (see lBlaizot et ai]l2006t) . However, a good agreement 
with the SPH simulation is achieved by simply associat- 
ing satellite galaxies with randomly selected subsets of dark 
matter halo particles or dark matter sub-structures. 



5 GALAXY GROWTH 

Galaxies grow by accreting gas and by merging with other 
galaxies. In this Section, we compare semi-analytic and SPH 
results at three different levels. First, we discuss the mass 
assembly of an individual galaxy (the central galaxy of the 
third most massive halo in the simulation at z = 0). We 
do this by following the "main branch" of the galaxy tree, 
both in the SPH and in the GallCS outputs. That is, we 
recursively link the galaxy to its most massive progenitor. 
Second, we compare semi-analytic and SPH predictions for 
the total cooling rate integrated over this halo's history. In 
the SPH simulation, we measure this quantity as the cu- 
mulative gas accretion rate of all the galaxy's resolved pro- 
genitors at a given time. In GallCS, we simply add up the 
mass of the gas that cools in all the progenitors of our halo 
at each time-step. Notice that these two estimates are dif- 
ferently affected by resolution effects. In the semi-analytic 
case, we only consider cooling in haloes that contain more 
than 20 dark matter particles. In the SPH case, the limit is 
more vague, because cooling is not affected by resolution in 
a straightforward way. Hence the comparison is not precise, 
but it suffices for qualitative assessment. Third, we compare 
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Figure 7. Gas accretion history of the galaxies in the third 
most massive halo contained in the computational box. The total 
amount of gas accreted by the galaxies in the halo at each cosmic 
time (black lines) has been decomposed into the contributions of 
the cold (blue) and the hot (red) mode. See Section 3.2 and Sec- 
tion 5 for an explanation of how we distinguish between "cold" 
and "hot" accretion. 



the total cooling rates in the entire simulation volume, for 
the SPH and the GallCS case. Here again, resolution affects 
the results slightly differently in the two methods. 

Figure 6 illustrates galaxy growth for the central object 
of the second most massive halo (Mh a i ~ 3 x 1O 13 M0). We 
start from the galaxy at z = and trace it back in time 
by selecting the most massive progenitor at each step. The 
black solid line shows the growth of its total galaxy baryonic 
mass in the SPH simulation. The no-feedback GallCS model, 
shown by the black dashed line, is in very good agreement 
with the SPH model. The only notable differences are due to 
merging. In the SPH model there is a sudden increase of the 
galaxy mass from 5 x 1O 1O M to 1.3 x 10 11 M at z ~ 3.3 - 
3.5, produced by a major merger between two outputs. The 
semi-analytic model contains a similar merging event but it 
occurs slightly earlier (between the snapshot at z ~ 3.7 and 
the snapshot at z ~ 3.5) and the growth is somewhat smaller 
(from 4 x 10 10 M Q to 8 x 10 10 Af Q ). The final difference in 
mass between the two galaxies is almost entirely due to a 
merger in the GallCS model at z ~ 0.039 - 0.076. In this 
merger, the twelfth since the first at z ~ 3.5 — 3.7, the galaxy 
mass has risen from ~ 7.4 x 10 n M« to ~ 9.5 x 10 U M B . 



GallCS 




Figure 8. The rate at which galaxies accrete gas in a cosmic 
comoving volume. The total value (black lines) is decomposed 
into the contributions of the cold mode (blue lines) and the hot 
mode (red lines). 



However, for the most part, the two growth curves nearly 
overlap. 

The red solid and dashed lines show the total amount of 
gas accreted by the galaxy since the start of the simulation 
in the two calculations. In the absence of mergers, the red 
and the black lines of each type would coincide. From the 
agreement between the red solid and dashed lines, we see 
that the two methods agree on the relative amounts of ac- 
cretion growth and merger growth. Finally, blue lines show 
the amount of cold gas in the galaxy as a function of red- 
shift. The two methods agree well down to z — 1.5. Below 
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this redshift, the SPH simulation predicts a lower gas frac- 
tion, but the gas fraction is small in each case. The residual 
gas fraction in this gas-poor regime depends on the details 
of the star formation prescription, so some disagreement is 
expected. 

iKatz et alJ <2003ft and lKeres et alJ <2005h have investi- 
gated gas accretion in the same SPH simulation that we are 
considering here. By tracing the temperature history of ac- 
creted particles, they have identified two distinct modes of 
gas accretion. About half of the gas follows the expected 
track in the conventional picture of galaxy formation: it 
shock-heats to the virial temperature of the galaxy's po- 
tential well (T ~ 10 6 K for a Milky Way type galaxy) before 
cooling, condensing, and forming stars. However, the other 
half radiates its acquired gravitational energy at much lower 
temperatures, typically T < 10 5 K. Cold accretion is often 
directed along filaments, allowing galaxies to draw gas ef- 
ficiently from large distances, and it is dominant for low 
mass galaxies (Mhalo < 3 x 1O 11 M0). Hot accretion is quasi- 
spherical and dominates the growth of high mass systems. 

iBirnboim fc Dekell fcOOfl and iDekel fc Birnboiml 

(2006) argue that the transition between the cold and hot 
accretion regimes is determined by the stability criterion 
for a shock at the halo's viral radius, which in turn depends 
on the ratio of the post-shock cooling tim escale to the 
dynam ical time. Motivated by this argument. ICroton et alJ 
(2005) suggest that the cold and hot accretion regimes in 
hydrodynamic simulations should be identified with the 
rapid cooling and slow cooling regimes in semi-analytic 
models. We investigate this identification in Figures 7 
and 8. 

The black curve in Figure 7a shows the total accretion 
rate onto all galaxies in the third most massive halo of the 
SPH simulation as a function of time. The blue and the 
red curves show the cold and hot accretion rates, respec- 
tively. Here cold accretion consists of gas particles whose 
maximum temperature at any phase of their evolution is 
Tmax < 2.5 x 10 5 K. B ecause the T max histogram is bi- 
modal llKeres et al .120051) . the division of hot and cold modes 
is insensitive to the exact choice of threshold. Figure 7b 
presents analogous results for the same halo in the no- 
feedback GallCS model. Here we identify cold accretion by 
the r coo i > rinfaii criterion discussed in Section 3.2. With this 
identification, the qualitative agreement between the SPH 
and GallCS calculations is quite good. In both cases there 
is a brief initial phase of cold accretion, when the universe 
is less than ~ 2 Gyr old, but hot accretion dominates after 
the halo becomes massive enough to support a virial shock. 
In the SPH simulation, cold and hot accr etion can co-exist 
for the same galaxy f see lKeres et alj|2005h . and in both cal- 
culations cold accretion continues onto lower mass galaxies 
that will subsequently join the main halo as satellites. How- 
ever, the central galaxy of this halo is massive (~ 1O 12 M0) 
and contains a large fraction of the cooled baryons, so it 
dominates the accretion statistics. The cold accretion rate 
on galaxies in this halo is negligible at t > 5 Gyr, in both 
calculations. 

Figures 8a and 8b present a similar comparison for the 
accretion rates averaged over the entire simulation volume. 
The results are analogous, but cold accretion remains impor- 
tant for a longer time because most haloes are less massive 



than the one shown in Figure 7. The agreement between the 
SPH and the GallCS calculations is now even better. 

Figure 8c shows the predictions of the full GallCS model 
with feedback. As discussed in Section 3.4, supernova feed- 
back is modelled by assuming that supernovae drive a wind 
with an outflow rate approximately equal to the star forma- 
tion rate, independent of the potential well's depth. Since 
ejected gas may cool and be accreted again, perhaps multi- 
ple times, the total accretion rate in the full GallCS model 
is higher than that in Figures 8a or 8b, even though the final 
galaxy masses are lower. Gas entrained in a galactic wind is 
heated to the halo virial temperature, but if its subsequent 
cooling time is shorter than the infall time, we count it as 
cold accretion, just as we do for infalling gas. As most super- 
nova feedback occurs in low mass galaxies with short cooling 
times, it is the cold accretion rate that rises significantly. 
Black hole feedback is modelled with a cooling cut-off in 
groups with EMbuige > 10 11 A/q, so it acts to suppress hot 
accretion in massive haloes. As a result, feedback reduces 
hot accretion in GallCS, especia l ly at late times. The new 
GallCS model of ICattaneo et ail l)200ft) makes the assump- 
tion that all hot accretion is suppressed by AGN feedback, 
though it incorporates a more sophisticated calculation of 
where the cold-to-hot transition occurs. 



6 THE EMERGENCE OF THE GALAXY 
BIMODALITY 

We now turn to our principal subject, the origin of bimodal- 
ity in the galaxy population. We must first understand the 
connection between galaxy gas content and star formation 
rate (which we denote either by SFR or by M*). Figure 9 
plots SFR against galaxy gas mass at z = 3 and z = 0, for 
the SPH, no-feedback GallCS, and full GallCS models. In 
each panel, red points represent satellite galaxies. Solid and 
dashed lines are the same in every panel and are provided 
for visual reference. 

At z — 3, the SPH and no-feedback GallCS populations 
show similar trends between gas mass and SFR, roughly 
SFR= 16(M gas /10 10 M Q ) 1 - 8 M Q yr" 1 . There is substantial 
scatter about the mean relation but no obvious difference 
between central and satellite galaxies of the same gas mass. 
The SFR at fixed gas mass in the full GallCS model is lower 
by a factor ~ 2, but the slope of the trend is similar. At 
z = 0, the SFR for a given M gas is lower than at z = 3, as ex- 
pected given the characteristically larger sizes and lower sur- 
face densities of low redshift galaxies. The SFRs are slightly 
higher at fixed M gas in the SPH simulation, but there are 
many fewer galaxies with large gas masses because massive 
galaxies in the SPH simulation convert nearly all of their gas 
into stars. The trends in Figure 9 are determined largely by 
the star formation law, which is qualitatively similar among 
the three models (even though the feedback prescriptions 
are different). 

Figure 10 plots gas fraction against galaxy baryonic 
mass. Here we see important differences among the three 
models and the first indications of bimodal galaxy popula- 
tions. As in Figure 9, blue points represent central galaxies 
and red points satellite galaxies. At z = 0, the no- feedback 
GallCS model shows two distinct, albeit fuzzy sequences in 
this Afgas/M ga i vs. A:f ga i plane. The gas rich sequence is cen- 
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Figure 9. The relation between gas content and star formation rate in the SPH simulation (top), in the GallCS model without supernova 
or AGN feedback (centre) and in the standard GallCS model (bottom). The blue points correspond to galaxies that are at the centre of 
their dark matter halo (isolated field galaxies, brightest group or cluster galaxies), while the red points are satellite galaxies in groups 
or clusters. The solid line and the dashed line correspond to l.SLogMg — 0.6 and 1.4LogMg — 1.1 , respectively (Mg = M gas /10 9 Mq). 
They are not fits and are just shown to guide the eye. 



tred on M g as/A/ ga i ~ 0.2, while the gas poor sequence is cen- 
tred on M gas /M ga i ~ 0.005, with a large scatter. Most cen- 
tral galaxies occupy the gas rich sequence and most satellite 
galaxies occupy the gas poor sequence, but there are excep- 
tions in both directions. The most obvious effect of feedback 
in the full GallCS model is to reduce the baryonic masses of 
the most massive galaxies, to a maximum of ~ 3 x 10 1 Mq. 
In contrast to the no-feedback GallCS case, many of the 
most massive central galaxies are gas poor instead of gas 
rich. 



In the SPH simulation, it is possible for a galaxy to 
have a gas fraction of zero, if all of the particles that com- 
pose it are purely stellar. We have arbitrarily assigned these 
galaxies a gas fraction of 10 -3 ' 5 with a 0.2 dex scatter, so 
that one can see them in Figure 10b. The dashed line in this 
panel shows the gas fraction for a galaxy that contains one 
gas SPH particle only. There is a significant gap between the 
bottom edge of the upper sequence of points and this dis- 
creteness limit. This gap suggests a genuinely bimodal be- 
haviour rather than an artificially divided continuum, with 
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Figure 10. The gas fraction as a function of galaxy mass and redshift. The blue points correspond to galaxies that are at the centre 
of their dark matter halo, while the red points are satellite galaxies. The diagonal line in (b) shows the position of galaxy with a given 
mass and one gas SPH particle only. The galaxies below this line do not contain any measurable mass of gas. They have been attributed 
a gas mass drawn from a 10 — 3 - 5 ± 2 log-normal distribution for the gas fraction to put them on the plot. 



all points above the line belonging to the gas rich sequence. 
Within this sequence, satellite galaxies have lower average 
gas fractions than central galaxies of the same mass. There is 
also a trend of decreasing gas fraction with increasing galaxy 
mass, but this is at least partly an artefact of resolution: gas 
densities are underestimated in low mass galaxies, and gas 
consumption by star formation is thus underestimated as 
well. Consequently, we do not regard this trend as a robust 
prediction, at least for masses below M ga i ~ 2 x 1O 1O M0 
(less than about 200 SPH particles). Once we account for 
the resolution and discreteness effects, the SPH results ap- 



pear qualitatively similar to the no-feedback GallCS results, 
with the notable exception that many SPH satellites occupy 
the low end of the gas rich sequence rather than the separate 
gas poor sequence. 

Figure 11 presents the principal result of this paper, 
the distribution of galaxies in the plane of stellar mass 
and star formation time-scale, M*/M*. Dividing the Hub- 
ble time to — 13.4 Gyr by this time-scale gives the ra- 
tio of the galaxy's current SFR to its time-averaged SFR. 
The light from galaxies with long star formation time-scales 
is dominated by red giants from old stellar populations, 
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Figure 11. The relation between galaxy mass and M*/M» at different redshifts in the SPH simulation (top), in the GallCS model 
without supernova or AGN feedback (centre) and in the standard GallCS model (bottom). M*/M t is the time that a galaxy would 
have taken to form its stars if it had formed them at its current star formation rate. A value of M„/M t of the order of the Hubble time 
suggests a constant star formation rate throughout the life of the Universe. A value much longer than the Hubble time implies that the 
star formation rate has been much higher in the past than it is now and that the current star formation is negligible. The blue points 
are galaxies at the centre of their halo and the red points are satellite galaxies. SPH galaxies with zero SFR have been displayed on the 
plot by attributing them a value of Log(M*/M*/Gyr) drawn from a 3 ± 0.2 Gaussian distribution. 



while young stars produce a significant amount of blue 
light in galaxies with more current star formation. Fig- 
ure 11 is thus a theoretical version of a galaxy "colour- 
magnitude" diagram. At low redshift, many SPH galaxies 
have zero star formation rate. We have arbitrarily assigned 
these galaxies time-scales of 10 3 Gyr, with random scatter 
of 0.2 dex, so that they appear on these plots. The Gal- 
lCS time-scale distributions have a sharp upper edge at 
M*/M* ~ 10 2 ' 4 Gyr. To understand its origin, consider the 
specific SFR, M,/M, = (M»/M gas )/(M gas /Af < ,). The first 
term is M,/M gas ~ 2(M gas /10 10 M Q )°- 4 M Q yr- 1 (Figure 9). 
The second term is M gas /M* > 10~ 3 in the GallCS model, 
but not in the SPH model, where many galaxies contain no 
gas at all (Figure 10). The tendency of GallCS galaxies to 
retain a minimum gas fraction explains the upper edge in 
the M*/M, distribution. 

The main features of Figure 11 follow naturally from the 
trends in Figures 9 and 10. At z = 0, the SPH simulation 



has a bimodal galaxy distribution. The old (red) sequence 
is populated mainly by satellite galaxies, with a few lower 
mass central galaxies mixed in. The young (blue) sequence 
contains a mix of central and satellite galaxies, but the con- 
sumption of gas in satellites is evidently fast enough to leave 
few galaxies with intermediate star formation time-scales. 
The most massive central galaxies all reside on the young 
sequence because of their continuing gas accretion, as seen 
already in Figure 10. Results for the no-feedback GallCS 
model are similar, except that there are fewer satellites in 
the young sequence. In GallCS, galaxies experience no gas 
accretion at all once they fall into larger haloes and become 
satellites. In the SPH simulation, some satellites continue 
to accrete gas, e.g. because they are the central objects of 
sub-haloes that have not yet dissolved and lost their iden- 
tity. The complete shutoff of accretion in GallCS satellites 
therefore appears overly drastic. 

Results for the full GallCS model are similar, with the 
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crucial distinction that many of the massive central galaxies 
are now on the old sequence. As we have already empha- 
sised, this observational success of the full GallCS model 
is a consequence of AGN feedback, which suppresses con- 
tinuing gas accretion in the most massive haloes. The SPH 
and no-feedback GallCS populations do include some fairly 
massive galaxies in the old sequence, but nearly all of these 
are satellites in the most massive halo (see Figure 2). This 
halo is anomalously large for this simulation volume, so a 
larger simulation would probably yield a smaller fraction of 
mas sive old galaxies in these two models. 

Kauffman n et alJ i2004T) have investigated the distribu- 
tion of M*/M* in the SDSS, using spectral line diagnostics 
to estimate star formation rates and stellar mass-to-light ra- 
tios. They obtain a bimodal distribution with peaks corre- 
sponding to log(M»/M»/Gyr) ~ 1.2 and 2.6, in reasonable 
agreement with the GallCS predictions in Figure 11. The 
typical star formation time-scales on the young sequence in 
the SPH simulation appear slightly too long. 

Moving to higher redshift, we see that all three models 
predict a unimodal distribution of star formation time-scales 
at z = 3, with typical value M*/M» ~ 1 Gyr. By z = 1, the 
time-scale on the young sequence has increased to about 
4 Gyr, and an old sequence has begun to emerge. In the SPH 
and no-feedback GallCS populations, this sequence mostly 
contains relatively low mass satellites. However, in the full 
GallCS model the central galaxies of the most massive z = 1 
haloes have already joined the "red and dead" population. 



7 SUMMARY AND DISCUSSION 

We have compared the galaxy po pulations of an SPH simu- 
lation, the hybrid GallCS model of lHatton et alJ <|2003T) . and 
a stripped down version of this model with no feedback from 
star formation or AGN. The GallCS calculations are applied 
to the dark matter halo population extracted from the SPH 
simulation, allowing direct comparison of the results. 

For the most part, the SPH and no-feedback GallCS 
predictions agree remarkably well. In particular, the two 
methods predict similar galaxy baryonic mass functions 
and similar dependence of these mass functions on envi- 
ronment and redshift. The SPH simulation predicts some- 
what more rapid growth at high redshift, probably a con- 
sequence of efficient filamentary accretion and a star for- 
mation time-scale that tracks the shorter dynamical time- 
scales of high redshift galaxies. The global gas accretion his- 
tories are similar in the two calculations, and this agree- 
ment extends to the individual galaxy level as shown in 
Figure 6. The breakdown between cold and hot accre- 
tion in the SPH simulation is qualitatively similar to the 
breakdown between the rapid cooling (r coo i > n n f a ii) and 
slow cooling (r coo i < r infa ii) accreti on in GallCS . This 
agreement supports the suggestion bv lCroton et alJ (l2005h 
that the cold /hot dichotomy seen in hydrodynamic simula- 
tions <Katz|l992l:lKatz et al.ll2003l: iBirnboim fc Dekel2003t 
iKeres et alJ l2005UDekel fc Birnboiml 12006ft can be identi- 
fied with the infall-dominated /cooli ng-dominated regimes 
of traditional semi-analytic models <W hite^^jVenkl f^9ll: 
iKauifmann et alll993tlCole et alJl994l) . lKeres et al.l (120051) 
note some important caveats to this identification, especially 



for models that incorporate a photoionising background, 
which increases cooling time-scales in low mass haloes. 

Both the SPH and the no-feedback GallCS models pre- 
dict a bimodal galaxy population with a "blue" sequence 
of gas-rich, star-forming galaxies and a "red" sequence of 
gas-poor galaxies dominated by old stellar populations. The 
red sequence is populated mainly by satellite galaxies, which 
stop accreting gas after they fall into larger haloes. A larger 
fraction of SPH satellites remain on the gas rich, star- 
forming sequence, which suggests that the truncation of gas 
accretion in GallCS is unrealistically abrupt. However, both 
calculations agree on a crucial point: the central galaxies of 
massive haloes experience continuing gas accretion, and they 
therefore remain on the star-forming sequence. The masses 
of these galaxies are substantially larger than those inferred 
for galaxies of similar space density in the real universe — 
i.e., both models drastically overpredict the high end of the 
galaxy baryonic mass function. There is also agreement that 
gas accretion by satellites in cluster-size haloes is negligible. 

While the agreement between SPH and no-feedback 
GallCS could in principle reflect a conspiracy of errors, it 
seems more reasonable to conclude that both methods do 
a fairly accurate job of modelling the growth of galaxies 
from cosmological initial conditions. This conclusion in turn 
implies that the observational failings of these models are a 
consequence of missing input physics. The full GallCS model 
is much more successful in reproducing observations because 
of its prescriptions for supernova and AGN feedback. Super- 
nova feedback reduces the masses of low and intermediate 
mass galaxies by about a factor of two, as the gas outflow 
rate is comparable to the star formation rate by construc- 
tion. AGN feedback produces a sharp cut-off in the bary- 
onic mass function by suppressing the cooling of gas in high 
mass haloes. Crucially, this mechanism also moves the cen- 
tral galaxies of massive haloes to the old, gas-poor sequence. 

Some mechanism that suppresses gas cooling in mas- 
sive haloes appears essential to explaining both the expo- 
nential cut-off in the observed galaxy lumin osity function 
jKauffmann et alJ Il993l: iBenson et all 120031) and the ob- 
served form of bimodality in the galaxy population, where 
nearly a ll of the most massive galaxies reside o n the red se- 
quence jBlanton et alJ l2003t iKauffmanrj 1200^ . A s argued 
bv iBirnboim fc Dekelf (I2003T) . iKatz et alJ (l2003ft. iBinnevI 
<2004h . IKeres et alJ J2005I) . and lDekel fc Birnboiml J2006|) . 
it seems natural to associate this mechanism with the "hot 
mode" of gas accretion that becomes dominant in haloes 
more massive than ~ 10 12 M Q , since the implied mass scale 
is close to t he tra nsition mass inferred observationally by 
Kauf fmannl i2003l) and since the amounts of cool gas in ob- 
served galaxy clusters are much lower than predicted by 
the X-ray cooling rates in the absence of heating (e.g., 
iKaastra et al .1200 ll) . In this view, the factor that determines 
the effectiveness of AGN feedback is not so much the mass 
of the black hole as the state of the gas that it will affect. 
Even a luminous AGN will have difficulty reversing the in- 
flow of dense, cold gas along filamentary streams, or ejecting 
the interstellar medium of a galaxy disc. However, an AGN 
can couple more easily to the diffuse, quasi-spherical hot at- 
mosphere of a massive halo, and, with plausible accretion 
rates and efficiencies, there is enough energy to suppress the 
gas coolin g that would oth erwise be expected in these haloes 
(see, e.g- lBest et"aT]l2006h . 
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In t his regard, it is i nteres ting that the new GallCS 
model of ICattaneo et all ]2006|) is con siderably more suc- 
cessful than the HHg^onelHiiT ||2P03) version of GallCS 
tested here. The Hatton et al. implementation ties AGN 
feedback to black hole mass and, thus, to bulge mass. This 
prescription suppresses cooling in massive haloes, but not 
effectively enough. Massive galaxies do lie on a "red" se- 
quence separated from the "blue" sequence of star-forming 
galaxies. Nevertheless, a quantitative comparison with the 
observations obtained using a stellar population synthesis 
model shows that they are still not red enough. Strong 
supernova feedback, achieved by assuming that the mass 
loss rate through a supernova-driven wind is approximately 
equal to the star formation rate independently of the depth 
of the potential well, can compensate the incomplete shut- 
down of cooling i n massive haloes, but it also ejects too 
many metals. The lCattaneo et al] (I2006T) model completely 
suppresses gas accretion in haloes above the critical mass 
where virial shocks aris e, computed via the methods of 
iDekel fc Birnboiml ]200d) . This implementation yields ex- 
cellent agreement with the observed colour-luminosity dis- 
tribut ion of galax i es at l ow and high redshift. 

ICroton et al] ]2005l) achieve a similar level of obser- 
vational success in a hybrid model that uses kinetic feed- 
back from radio-loud AGN to suppress gas cooling in mas- 
sive haloes. Because this mechanism is assumed to affect 
only the quasi-hydrostatic atmospheres that arise in haloes 
with long post-shock cooling times, the practical effect of 
this scheme is quite similar to the sharp c ut-off in the 
ICattaneo et al.l i2006l) model. Very recently, iBower et al] 
(2006) ha ve presented an o ther m odel that is very similar 
to that of ICattaneo et al] M2006T) . both in the method and 
the key results. 

There is still much to be understood about the physics 
that suppresses gas cooling in massive haloes, i.e., the de- 
tailed mechanism by which AGN energy couples to the sur- 
rounding gas, or even whether AGN are really the source 
of this suppression at all. However, it is striking that the 
theoretically derived mass scale separating the cold and hot 
accretion regimes corresponds so well to the observed mass 
scale that marks the transition between the blue and red 
galaxy populations, and that the most successful models of 
the observed galaxy population are those that incorporate 
suppressed hot accretion to identify one transition with the 
other. 
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